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Abstract 

We investigate the quantum phase transition in the random transverse-field Ising model under the influence of Ohmic 
dissipation. To this end, we numerically implement a strong-disorder renormalization-group scheme. We find that Ohmic 
dissipation destroys the quantum critical point and the associated quantum Griffiths phase by smearing. Our results 
quantitatively confirm a recent theory [Phys. Rev. Lett. 100, 240601 (2008)] of smeared quantum phase transitions. 
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1. Introduction 

The presence of impurities, defects, or other types of 
quenched disorder can qualitatively change the proper- 
ties of a continuous phase transition. At thermal tran- 
sitions, this is controlled by the Harris criterion [1]: If 
a clean critical point violates the inequality v > 2/d 
where v is the correlation length exponent and d is 
the space dimensionality, it is perturbatively unstable 
against weak disorder. The generic result of adding dis- 
order to such a system is a new critical point with dif- 
ferent critical exponents but conventional power-law 
scaling. Moreover, Griffiths [2] showed that the free en- 
ergy is singular not just at the transition, but in an en- 
tire parameter region around the transition now known 
as the Griffiths phase. This is caused by rare spatial 
regions that are locally in the ordered phase while the 
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bulk system is still in the disordered phase. However, it 
was soon realized that thermodynamic Griffiths effects 
are generically very weak and probably unobservable 
in experiment [3]. 

At zero-temperature quantum phase transitions, 
disorder effects can be dramatically stronger than 
at thermal transition. Fisher [4,5] showed that the 
random transverse-field Ising chain has an exotic 
infinite-randomness critical point featuring ultraslow 
activated rather than power-law dynamical scaling. 
The associated quantum Griffiths singularities take 
power-law forms [6,7], implying a diverging suscepti- 
bility in the Griffiths phase (for a review on rare region 
and Griffiths phenomena see, e.g., Ref. [8]). 

Dissipation can further enhance the effects of dis- 
order on a quantum phase transition. For Ising order 
parameter symmetry, each locally ordered rare region 
acts as two-level system which can undergo the local- 
ization transition of the spin-boson problem when cou- 
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pled to an Ohmic bath [9]. Thus, each region freezes 
independently of the bulk system, destroying the sharp 
phase transition by smearing [10]. 

Recently, we developed an analytical strong-disorder 
renormalization-group theory for the dissipative ran- 
dom transverse-field Ising chain that allowed us to de- 
termine the low-energy fixed points exactly [11]. We 
found that Fisher's infinite randomness critical point 
[4,5] and the quantum Griffiths phase are indeed de- 
stroyed by the dissipation. Instead, there is only one 
nontrivial fixed point describing an inhomogeneously 
ordered phase (the tail of the smeared transition). 

In this paper, we report the results of a numerical 
implementation of the strong-disorder renormalization 
group for the dissipative random transverse-field Ising 
chain. The purpose of the work is twofold; (i) it allows 
us to confirm and illustrate the theoretical predictions 
for the asymptotic low-energy behavior, (ii) More im- 
portantly, the numerical simulations allow us to test 
whether moderately disordered realistic systems actu- 
ally flow to the predicted fixed points, and they allow 
us to find the relevant crossover scales. The paper is 
organized as follows. We introduce the model Hamil- 
tonian in Sec. 2. We summarize the main results of the 
analytical renormalization-group theory in Sec. 3. Sec- 
tion 4 is devoted to the computer simulations. Finally, 
we conclude in Sec. 5. 



2. Dissipative random transverse-field Ising 
model 

We consider a one-dimensional random transverse- 
field Ising model, a prototypical model displaying an 
infinite-randomness quantum critical point. Dissipa- 
tion is introduced by coupling each spin linearly to an 
independent Ohmic bath of harmonic oscillators. The 
resulting Hamiltonian reads 

h = — Jiai<Ji +1 - ^2 

i i 
i,n i,n 

Here, a^' z are the Pauli matrices representing the spin- 
1/2 at site i. The bonds J; and transverse fields hi 
are independent random variables with bare (initial) 
probability distributions Pi (J) and Ri(h). a\ n (aj,„) 



are the creation (annihilation) operators of the n-th 
oscillator coupled to spin Oi via Aj, n , and Vi^ n is its 
frequency. All baths have the same bare Ohmic spectral 
function 

£(u) = 7r^ \% n 5{uj - Vi, n ) = 2 7 raa;e"" / " < = , (2) 

n 

with a the dimensionless dissipation strength and cj c 
the (bare) cutoff energy. Under the renormalization 
group, the cutoff will change, and the dissipation 
strength will become site-dependent. 



3. Strong-disorder renormalization group 

In this section, we briefly summarize our analytical 
strong-disorder renormalization-group approach [11] 
to the Hamiltonian (1). It is related to a numerical 
scheme suggested by Schehr and Rieger [12,13]. How- 
ever, treating the oscillator modes on equal footing 
with the spin degrees of freedom allows us to solve the 
problem analytically. 

The basic idea of the strong-disorder renormaliza- 
tion group [14,15] is to successively integrate out local 
high-energy states. In the Hamiltonian (1), the relevant 
local energies are the transverse fields hi, the bonds Ji, 
and the oscillator frequencies u i%n - In more detail, the 
renormalization group proceeds as follows. We start by 
determining the energy cutoff, i.e., the largest local en- 
ergy in the system Q, — max(/ij, Ji,iv c /p) where p>l 
is an arbitrary constant. In each renormalization-group 
step we now reduce the energy scale from $1 to — dQ 
by (i) integrating out all oscillators (at all sites i) with 
frequencies between p(fl — dfl) and pQ, using adiabatic 
renormalization [9] and (ii) decimating all transverse 
fields and all bonds between (Q — dfl) and in pertur- 
bation theory. 

In step (i) all transverse fields renormalize according 

to 

hi = hi (l - a i-^Pj (3) 

while the bonds remain unchanged. Here a; is the 
renormalized dissipation strength at site i. In step (ii), 
if two sites are coupled by a strong bond Ji — £1, the 
spins Oi and a i+ i can be treated as an effective spin 
cluster a with moment 

/i = m + fii+i , (4) 
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in a bath of renormalized dissipation strength 

a = on + a»+i = a(fn + fn+i) = ap, , (5) 

and renormalized transverse field 

h = hih i+ i/Ji . (6) 

Conversely, if a site experiences a strong field hi = fl, 
the corresponding spin ai is eliminated, creating a new 
bond between sites i — 1 and i + 

J — J i ^ 1 J i /h i . (7) 

We now iterate the complete renormalization-group 
step (3)-(7) decreasing the cutoff energy scale Q. Us- 
ing logarithmic variables V = In(f2j/f2) [where Qi 
is the initial (bare) value of fl], £ = ln(fi/J) and 
/3 = ln(fi//i), we can derive renormalization-group flow 
equations for the probability distribution V(Q of the 
bonds and the joint distribution 1Z(/3, /i) of the fields 
and moments. They are given by 

dV dV ( C \ 

If = 7^ + (i-«fc>)^ (o) (p®p) 

+ [P(0)-(l-a/Z )7^(0)]P, (8) 

f = d-^)f + 7>(0)(tt#tt) 

-[7>(0)-(l-a/Z )7^(0)]ft , (9) 

where 1Z p (/3) = J °° /i)d/i is the distribution of the 
fields and jS is the average moment of clusters about to 
be decimated (defined by ~p 7Z.p (0) = / °° /j,TZ(0, /J,)d[i). 

The symbol V ®V = ^V{C!)V{C, - C'X denotes 
the convolution. The first term on the r.h.s. of (8) and 
(9) is due to the rescaling of £ and j3 with T and the 
renormalization (3) of h by the baths. The second term 
implements the recursion relations (4), (6) and (7) for 
the moments, fields and bonds. The last term ensures 
the normalization of V and 1Z. As expected, for a = 0, 
(8) and (9) become identical to the dissipationless case 
[4,5]. 

The qualitative change of the physics brought about 
by the dissipation can be seen already in the probability 
of decimating a field, (f — afI )TZp(0), which decreases 
with increasing dissipation strength and cluster size. 
Clusters with moment [x > 1/a are never decimated. 
Thus, in the presence of dissipation, the flow equations 
contain a finite length scale above which the cluster 
dynamics freezes. 



In Ref. [11], we determined the fixed points, i.e., sta- 
tionary solutions of the flow equations (8) and (9) (in- 
variant under a general rescaling r\ — C//c(r), = 
/3/fp(T) and v = /u// p (r)). They correspond to stable 
phases or critical points. We found that the dissipation 
destroys the the critical fixed point of the dissipation- 
less model [4,5] and the ones associated with the cor- 
responding disordered and ordered quantum Griffiths 
phases. Instead, for overlapping bond and field distri- 
butions, there is only one line of well-behaved fixed 
points (parameterized by Va > 0) corresponding to the 
tail of the ordered phase. Here, /<; = 1, = exp('Por), 
fp — Texp('Por)- The fixed-point distributions can be 
found in closed form, they read 

V*(0=Voe- Vo< , (10) 
TZ*(e, v) = TZae-^SiO - av) , (11) 
with 72-0 being a nonuniversal constant. This fixed point 
is similar to the ordered Griffiths phase in the dissipa- 
tionless case, but fp // M — > oo as r — > oo. Transforming 
the field distribution (11) back to the original trans- 
verse fields h gives power-law behavior ~ h 710 ^ 01 ^^ 1 . 
This result implies that there is no fixed-point solution 
with /<; / fp — > const as Y — > oo in the presence of dissi- 
pation, proving that there is no quantum critical point 
where fields and bonds compete at all energy scales. 

4. Numerical results 

In this main section of the paper, we present nu- 
merical renormalization-group results for the Hamilto- 
nian (1) with moderate disorder. The goal is to check 
whether the renormalization-group flow is indeed to- 
wards the fixed point (10,11) and to investigate the 
crossover from dissipationless to dissipative behavior 
for small dissipation strength. 

To do so, we numerically implement a strong- 
disorder renormalization-group scheme similar to 
Refs. [12,13], and apply it to very long chains up to 
Lo = 10 6 sites. All data are averages over 10 3 different 
disorder realizations. The random bonds and fields are 
drawn from probability distributions 

Pl(J)xJ- X (Jmin < J < Jmax) , (12) 

Ri(h) oc h~ x (h min <h<h 

max 

) • (13) 

All sites have the same initial magnetic moment fn = 1 
and oscillator cutoff uj c = pftmax. We have used p = 1, 
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Fig. 1. Renormalization-group flow of (f) (open symbols) 
and 5£ (filled symbols) as a function of the cutoff energy 
scale, r = In(f2//Q) for several different chains. Note that 
(f) and <5£ agree to good approximation. 

but because p does not appear in the recursion rela- 
tions (3)-(7), using any other value just amounts to a 
redefinition of the bare distribution Ri(h). Using this 
method, we have investigated a large number of differ- 
ent parameter sets. In following we present a selection 
of typical results. 

Let us start by considering the renormalization- 
group flow of the averages and standard deviations 
of the logarithmic bond and field variables as well as 
the magnetic moment. According to the fixed-point 
solution (10,11), they are expected to behave as 



(C) = 5( = 1/Vo , 
(/3) = 5/3 = aTe v ° T /1l , 
(jj) =8^ = e v ° r /TZ 



(14) 

(15) 
(16) 

in the low-energy limit F — > oo. In Fig. I, we show 
the evolution of the average and standard deviation of 
the logarithmic bond variable £ for several chains with 
different /i max and J max - 2 The dissipation strength 
is fixed at a = 0.1. Without dissipation, the chain 
with (ft-max, Jmax) = (0.9, 1.0) would be in the ordered 
phase, the chain with (1.0,1.0) would be critical, and 
all others would be in the disordered phase. The figure 
shows that (£) and <5£ initially increase under renor- 
malization (corresponding to a rapid drop of the bond 
energies J); however after a sharp crossover, they set- 
tle on a constant value as expected from (14). 

Figure 2 shows the corresponding data for the av- 
erages of the logarithmic field variable f3, the cluster 



2 Before measuring any quantity we always integrate out 
the oscillators such that all local cutoffs ui c i are identical. 



Fig. 2. Renormalization-group flow of (/3), (fi) (open sym- 
bols) and <5/3, 8/i (filled symbols) as a function of the cut- 
off energy scale, V = In(f2//f2). The solid lines represent 
and 5/3 while the dashed lines represent (fi) and <5/i. 
Again, averages and standard deviations agree very well for 
T -> oo. 

magnetic moment /i, and their respective standard 
deviations. They are plotted as \n({/3) / (aF)) , ln({/i)), 
ln(<5/3 / (ar)), and ln(<5/i) versus T. In this plot, the 
fixed-point forms (15), (16) correspond to straight 
lines. The figure shows that all data indeed follow the 
expected behavior after some short initial transients. 

So far, we have analyzed the renormalization-group 
flow of the averages and standard deviations of the 
bond and field variables as well as the magnetic mo- 
ments. For a full confirmation of the theoretical predic- 
tions (10) and (11), we need to analyze the full prob- 
ability distributions. Figure 3 shows snapshots of the 
probability distributions P(C) ; TZe(@) and TZ u (v) along 
the renormalization-group flow. Here, we have used the 
rescaled field variable 9 = j3/{Te p ° T ) and the rescaled 
moment v = fj,/(re p °). The figure shows that all dis- 
tributions quickly approach the predicted functional 
forms. In the late renormalization-group stages, the 
numerical data are in excellent quantitative agreement 
with the theoretical results (10) and (11). 

After having analyzed the behavior at fixed value 
of the dissipation strength a, we now turn to the a- 
dependence of the renormalization-group flow. Figure 
4 shows the evolution of {() , (/3), and (fi) as well as their 
standard deviations for different a. Because /i max = 
1.0 is larger than J ma x = 0.8, the dissipationless chain 
(q = 0) is in the disordered (Griffiths) phase. This can 
be seen from the rapid increase of the bond variable 
((), corresponding to a rapid drop of the interaction 
energies J under renormalization while the field vari- 
able (P) approaches a constant and the moment (/z) of 
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Fig. 3. Probability distributions V{Q, Kg{0) and H v [y 
taken at different stages of the renormalization-group flow 
parametrized by the average distance L of the remaining 
spin clusters. The dashed lines represent the fixed-point 
distributions (10) and (11). (/i max , Jmax) = (1.0,0.8), the 
parameters are as in Fig. 1 and 2. 

the surviving clusters increases linearly with F [4,5]. 

In the presence of dissipation, the flow changes qual- 
itatively when the typical cluster size reaches 1/q, as 
can be seen in the third panel of Fig. 4. Beyond this 
crossover scale, it is now the bond variable {() that 
saturates while the field variable (/3) and the moment 
{(j,) rapidly increase, in agreement with Eqs. (14), (15), 
and (16). Thus, at a cluster moment of 1/a the flow 



Fig. 4. Renormalization-group flow of (£), ((3), and (fi) 
(open symbols) as well as their standard deviations 
(filled symbols) as a function of the cutoff energy scale, 
T = ln(Q//n) for (h Jmax) = (1.0,0.8). 

character changes from that of the disordered Griffiths 
phase to that of our inhomogeneously ordered phase. 
This crossover is caused by the fact that clusters with 
ctfj, > 1 undergo the localization transition of the spin- 
boson problem [9] and are not decimated under further 
renormalization. 

Finally, we have also studied the ultimate fate of the 
pseudo-critical point identified at intermediate scales 
by Schehr and Rieger. To this end, we have repeated the 



5 



< 
cf 

a 



-10, 







a=0.052 
/! max =l 




-L=f 






7 max =0.34 




-L=2 1 
~L=2 S 






m/p=10 4 




L=2 9 








- L=2 10 
L=2" 




/ // \ \ \ 






- L=2 12 
L=2 13 
-L=2 14 













*J -6 



10 



20 30 40 
InCfij/A) 



50 



60 



-10, 





-L=2 5 




-L=2 6 




-L=2 7 
-L=2 8 
L=2' 




- L=2'° 


' J a=0.052 


L=2" 


-J h im =i >^r-^^ 


- L=2' 2 


■ I ■ / m ! « =0 ' 34 ^VT^~\ 


L=2 13 


-1 ro/p=l0 4 NX 


-i=2" 




l|/=0.32 \ 





10 



ln(iyA)/L T 



Fig. 5. Probability distribution of the smallest energy gap A 
of a finite-size chain for different chain lengths. The param- 
eters /l max = 1.0, Jmax = 0.34, /l mill = J min = 0, X = 0, 

and a = 0.052 correspond to the pseudo-critical point of 
Ref. [12]. 

calculation of Ref. [12] for much longer chains of up to 
16000 sites (averaging over 10 6 disorder realizations). 
The quantity studied is the smallest excitation energy 
A of a finite chain, as estimated by the last nonzero 
transverse-field in the renormalization procedure. The 
upper panel of Fig. 5 shows the probability distribution 
of In(f2i/A) for different chain sizes. Following Ref. 
[12], only chains that are not yet frozen (A ^ 0) at the 
end of the renormalization procedure are included in 
the distribution. The parameters ft max = 1.0, J max = 
0.34, x = 0, and a — 0.052 exactly correspond to the 
pseudo-critical point of [12]. The figure shows that the 
character of the distribution changes significantly with 
increasing chain length L. The lower panel demon- 
strates that the data for L J$ 512 and not too large 
ln(fij /A) can be approximately scaled according to the 
pseudo-critical activated scaling form ln(S7r/A) ~ V* 
with ip — 0.32. However, for longer chains the distribu- 
tions do not scale. Instead, they become much broader, 
indicating a crossover to the functional form (11) char- 



acterizing the inhomogeneously ordered phase. Thus, 
the pseudo-critical behavior applies only to a transient 
regime of the renormalization-group flow and thus only 
to a transient energy or temperature window [13]. 



5. Conclusions 

In summary, we have presented results of an exten- 
sive numerical strong-disorder renormalization group 
for the random transverse- field Ising model with Ohmic 
dissipation. Our simulations quantitatively confirm the 
analytical theory of Ref. [11]. Specifically, we have ver- 
ified that the Ohmic dissipation destabilizes the quan- 
tum critical of the dissipationless system and the asso- 
ciated quantum Griffiths phase. No other critical point 
has been found. Instead, the low-energy behavior in 
the region of overlapping field and bond distributions 
is governed by a new line of fixed points describing 
the inhomogeneous tail of the ordered phase. Thus, the 
sharp quantum phase transition is destroyed by smear- 
ing due to the interplay of disorder and Ohmic dis- 
sipation. Note that Ohmic dissipation also suppresses 
the quantum Griffiths singularities at the percolation 
quantum phase transition [16] in a diluted transverse- 
field Ising model [17]. However, the percolation transi- 
tion remains sharp because it is driven by the critical 
geometry of the lattice. 

In addition to confirming the fixed-point structure 
of the analytical theory [11], our extensive numerical 
results also show that moderately disordered systems 
generically flow towards the new line of fixed points. 
The crossover to the asymptotic behavior occurs when 
the typical cluster moment reaches 1/a. 

Let us conclude by putting our results into a broader 
perspective. Recently, a general classification has been 
put forward of phase transitions in the presence of weak 
disorder [8] . It is based on the effective dimensionality 
of the defects or rare regions. Three classes need to 
be distinguished, (i) If the defect dimension is below 
the lower critical dimension dJT of the problem, the 
behavior is conventional; (ii) if it is right at , the 
transition is of infinite-randomness type; and (iii) if 
it is above d~ , finite clusters can order independently 
leading to a smeared transition. In our case, individual 
rare regions can undergo the localization transition of 
the spin-boson problem [9]. The system therefore falls 
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into the smeared-transition class. 

The results for a dissipative Ising magnet must be 
contrasted with the behavior of systems with continu- 
ous O(N) symmetry. While large Ising clusters freeze 
in the presence of Ohmic dissipation, O(N) clusters 
continue to fluctuate with a rate exponentially small 
in their moment [18], putting the system into class 
(ii). This leads to a sharp transition controlled by an 
infinite-randomness critical point in the same univer- 
sality class as the dissipationless random transverse- 
field Ising model [19,20]. 

Our results directly apply to quantum phase tran- 
sitions in disordered systems with discrete order 
parameter symmetry and Ohmic dissipation. The 
renormalization-group approach should be broadly 
applicable to a variety of disordered dissipative quan- 
tum systems such as arrays of resistively shunted 
Josephson junctions. 
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